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Abstract 

Two-pion induced self-energy contributions to the /O-meson mass are exam- 
ined in relation to the quenched approximation of QCD, where the physics 
associated with two-pion intermediate states has been excluded from vector- 
isovector correlation functions. Corrections to quenched QCD calculations of 
the yO-meson mass are estimated to be small at the order of a few percent of 
the /9-meson mass. The two-pion contributions display nonanalytic behavior 
as a function of the pion mass as the two-pion cut is encountered. The impli- 
cations of this nonanalytic behavior in extrapolations of full QCD calculations 
are also discussed. We note that for full QCD, the error made in making a 
linear extrapolation of the p mass, neglecting nonanalytic behavior, increases 

as one approaches the two-pion cut. 
12.38.-t, 12.38.Gc, 14.40.-n, 14.40.Cs 
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I. INTRODUCTION 



The lattice regularized approach to quantum field theory provides the best forum for the 
examination of the fundamental nonperturbative aspects of QCD. In the low momentum 
transfer regime, it is the only approach which in the foreseeable future holds a reasonable 
promise of confirming or rejecting the validity of QCD as the underlying theory of the strong 
interactions. 

Recently, it has become possible to perform quenched QCD calculations in which all 
systematic uncertainties are quantitatively estimated. Thus if the effects of quenching can 
be understood, the validity of QCD may be tested in the nonperturbative regime. Of 
particular note is the recent determination of the QCD coupling constant, otjjs^ from the 
IS — IP mass splitting of charmonium In this case one believes that the effects of 
quenching may be estimated with minimal model dependence. Corrections to the utilization 
of the quenched approximation of QCD have been estimated and are currently the dominant 
source of uncertainty in the final predictions 0. 

In this paper we will continue efforts along this line through the examination of systematic 
uncertainties in hadron mass spectrum calculations. In particular, the importance of the 
two-pion induced self-energy contribution to the p-meson mass is evaluated in relation to 
the quenched approximation of QCD and to full QCD. In the quenched approximation, the 
physics associated with two-pion intermediate states has been excluded in the numerical 
simulations. 

This investigation is motivated by recent results from the GFll group for the low-lying 
hadron mass spectrum in the quenched approximation of QCD. Their analysis is the first 
to systematically extrapolate QCD calculations to physical quark mass, zero lattice spacing 
and infinite volume. Their predictions display an impressive agreement with experiment. 
Of eight hadron mass ratios, six agree within one standard deviation and the remaining two 
ratios agree within I.60". However, since the quenched approximation leaves out so much 
important physics, one might question whether these results are actually too good [^] . 
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In the perturbative regime, many of the effects of not including disconnected quark loops 
when preparing an ensemble of gauge configurations may be accounted for in a simple renor- 
malization of the strong coupling constant. However, one also anticipates nonperturbative 
effects in making the quenched approximation. Unlike a global renormalization of the cou- 
pling constant, these effects are expected to be channel specific. For example, the quenched 
approximation of QCD leaves out the physics associated with the decay of the p meson to 
two pions. This physics must be accounted for and added to the quenched results prior to 
comparing with experimental data. Moreover, the calculated hadron masses are extrapo- 
lated as a function of the pion mass squared to the point at which the pion mass vanishes 
using linear extrapolation functions. Such an approach neglects nonlinear and indeed, non- 
analytic behavior in the continuum extrapolation function. For example, in the case of the 
p- meson mass, one expects nonanalyticity associated with the onset of the two-pion cut. 

A priori one does not know the relative importance of two-pion intermediate states of 
the p meson in describing the p mass. The substantial width of the p meson at 151.5 ± 1.2 
MeV indicates its coupling to pions is not small and correspondingly these dynamics may 
have significant influence on the p-meson mass. 

Geiger and Isgur were the first to study the possible importance of the two-pion induced 
self-energy of the p meson in relation to lattice QCD calculations P|. Their results are based 
on a string breaking quark model and predict large corrections to smooth extrapolations of 
the p-meson mass at approximately 70 MeV. These authors were motivated by the long 
standing problem of QCD predictions for the A^/p-mass ratio being too large. Their hope 
was the two-pion induced self-energy correction would sufficiently raise the p mass to solve 
this problem. However, we now understand that both the finite lattice spacing and the finite 
volume of the lattice act together to push up this mass ratio. Current estimates [Q] for this 
ratio corrected to the infinite volume, continuum limit are 1.22 ±0.11 in excellent agreement 
with the experimental value of 1.22. 

Geiger and Isgur |]^ advocate using the nonlinearities in the p mass as a function of the 
quark mass to correct for the linear extrapolation of lattice results. We wish to stress that 
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such a procedure is only sensible for the extrapolation of full QCD calculations. As we shall 
argue, the entire two-pion induced self-energy is absent in the quenched approximation and 
should be added on to quenched QCD results prior to comparison with experiment. We 
note that for the model of Ref . , quenched lattice calculations of the p-meson mass would 
be reduced by 160 MeV instead of increased by 70 MeV. This would further exacerbate the 
^^iriM/mp problem" discussed in their paper, rather that curing it. 

In the quenched approximation, the p meson cannot decay to two light pseudoscalars. 
As discussed in Ref. it is not possible to generate intermediate states of the p meson in 
the quenched approximation in which one has two isovector pseudoscalars. One might worry 
about the presence of the light isoscalar pseudoscalar rj' which fails to obtain its heavy mass 
in the quenched approximation [§-0. However, decay of the p meson to a vr 77' is forbidden 
by G-parity and decay to two ?7'-mesons is of course forbidden by the isospin invariance of 
the strong interactions. Hence the physics associated with two-pion intermediate states of 
the p should simply be added onto the results extracted from calculations of quenched QCD, 
provided the lattice spacing is defined by physical observables which do not have a similar 
dependence on pion decays. 

Current calculations of full QCD typically employ quark masses which place 2m^ > 
TUp. As a result, the functional form of the extrapolation function should account for the 
nonanalytic behavior in the p-meson mass as the two-pion cut is encountered. On the lattice, 
the spectral density does not have a cut but rather a series of poles at the points satisfying 

^s = 2{ml+pl)"\ (1) 

where pn are the discrete momenta allowed on the lattice. Obviously, to fully account for 
the two-pion induced self-energy of the p, one must first extrapolate the lattice results to 
zero lattice spacing and infinite volume prior to extrapolating the quark masses to physical 
values. 

In evaluating the integrals describing the coupling of pions to the p one must take into 
account the dependence of the p to two vr coupling constant Qp.,,.,, refiecting the internal 
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structure of these mesons. While this dependence was extracted from a string breaking 
quark model in Ref. 0, we elect to take a more agnostic approach and consider different 
methods of cutting off the integral. In particular we investigate a sharp ^-function cutoff and 
a dipole cutoff. The underlying reason for our agnosticism is the belief that models of the 
sort underlying Ref. ^ are not likely to correctly describe the structure of pseudo-Goldstone 
bosons such as the pion, as they do not incorporate chiral symmetry. While there are models 



of the pTTTT vertex which incorporate chiral symmetry and chiral symmetry breaking |[1(]|-[12 



these models are based on particular dynamical assumptions. Accordingly, it is difficult to 
assess the reliability of such models. Instead, we consider a range of possibilities for the 
vertex and to simplify this task, we consider convenient phenomenological forms. 

One other paper addressing this issue |jl3| sidesteps the problems surrounding the de- 
pendence of Qp-^T, by fixing gp^^T, to a constant and making two subtractions of the divergent 
integral at = 0. These subtractions are absorbed into a mass and wave function renormal- 
ization. However, this approach excludes any analysis of p-meson mass extrapolations as the 
subtraction terms themselves have an unknown dependence which has been lost in the 
renormalization procedure. Moreover, contributions from virtual two-pion states have been 
absorbed into the bare lattice parameters which is inconsistent with the dynamics contained 
in the quenched approximation. 

The outline of this paper is as follows. In Section II the model used in examining the 
two-pion induced self-energy is outlined. Two methods for regulating the divergent self- 
energy are explored. In section III the relevance of the self-energy corrections to quenched 
QCD simulations is discussed. Section IV addresses the quark mass extrapolation of full 
QCD calculations and the importance of nonlinear behavior in the p-meson mass. Finally, 
the implications of this investigation are summarized in Section V. 
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II. THE SELF-ENERGY 



In modeling the two-pion induced self-energy of the p meson, Sp^^, the standard pirn 
interaction motivated by low-energy current algebra is used. The effective Lagrange inter- 
action has the form [jl4| 

>Cint = -i gpTvn [t^ 7?j + gj^^ TT^ . (2) 

The pions are further assumed to interact exclusively through the p channel as summarized 
in the following Schwinger-Dyson equation for the p propagator 



(3) 



where 
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in Landau gauge, and Mq is the bare p-meson mass. The self-energy Sp^^ is defined through 
the solution 
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where 
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S*^"^ is given by the standard one loop integrals 
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Physically, the integral is convergent due to the momentum dependence of gpmr- However, 
this momentum dependence is unknown. In light of this uncertainty, it is reasonable to 
parameterize the momentum dependence in terms of some given functional form with an 
adjustable parameter controlling how the function falls off as a function of momentum trans- 
fer. We consider two regulation prescriptions. In one, we assume a monopole form for gpnn, 
and for comparison, we also consider a sharp 6'-function cutoff. 
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The simplest fashion for introducing a covariant cutoff function is through the use of a 
dispersion relation. The second term of (|^) is q independent and serves only to subtract the 
quadratic divergence of the first term in maintaining current conservation. As a result, we 
write a dispersion relation for S(g^) with one subtraction at = 0, 

s(,^),ir<i.,?!i!!iE£). (8) 

TT JO s s — 

Of course, the imaginary part of SpTrvr may be easily determined using any number of tech- 
niques for rendering the integral of finite. The imaginary part is 

= |j (l - (q- -Ami). (9) 

The value of gp.„T, at q^ = rrij^ is fixed by equating the imaginary parts of 



+ Ep^^ =[mp + — ) , (10) 



at q = m . The physical values [T^ nip = 768.1 MeV and T = 151.5 MeV fix gp^.^^ at ~ 6.0. 



A. ^-Function CutofF 

To illustrate the physics associated with the real part of the self-energy we first consider 
the integral of (H) cut off covariantly by a sharp ^-function at s = A^. The functional form 
is 



where 



.,^(l-l^)"\and..^(l-l^)"\ (12) 

Figure |1| illustrates the real part of Sp^^r evaluated at = for a variety of cutoffs 
ranging from slightly above to 4 GeV^. For small cutoffs, most of the strength in the 
integral lies below the p mass and consequently the p mass is pushed up due to the mixing 
with pion states. Of course, this behavior is completely consistent with that anticipated by 
simple quantum mechanical arguments. For larger cutoffs the strength above the p mass 
acts to reduce the p mass. 



B. Dipole Cutoff 



While the ^-function is useful as an illustrative tool, it suffers from being physically 
artificial and the results can be very sensitive to the value of the cutoff, as illustrated in 
figure |T[ In an attempt to better represent the dependence of (^ptttt a monopole form for 



This form maintains the normalization of (^ptttt defined at the physical p mass and renders 
the integral of (|[) finite. 

Of course, this approach is not without a few unphysical side effects. The most obvious 
problem is the introduction of spurious poles in the space like region for the s dependence 
of (y'pTTTr- However, the dispersion integral only samples the time like region and the presence 
of these unphysical poles should not affect the results. One could consider other functional 
forms. However, the effective physical value for the regulator mass. A, is itself unknown. 
Our aim is to estimate the importance of the two-pion induced self-energy relative to the p 
mass, as opposed to attempting to evaluate the actual correction. For this reason we view 
a consideration of the dipole regulator to be adequate. 

Evaluation of the dispersion relation of with (|T3|) leads to the following functional 
form for the real part of Sp^rTr 



each vertex is introduced and the dispersion relation of is evaluated with 




(13) 




(14) 



where 




(15) 



The imaginary part is recovered as in (|^). 
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Figure g illustrates the real part of the self-energy and its derivative with respect to 
at = 2 GeV^. The derivative clearly displays the nonanalytic behavior encountered at 
rrip = 2m7r. The second derivative is discontinuous at rrip = 2m7r and is infinite from above. 
The imaginary part of the self-energy is also illustrated in figure 0. 

A comparison with figure |I| indicates that at A^ = 2 GeV^ the results are not too sensitive 
to the manner in which the integral is regulated. Figure ^ illustrates the real part of the 
self-energy for the same values of A^ used in figure |l|. The sensitivity of the results to the 
value of A is greatly reduced and all curves display the same qualitative behavior. 

In the limit of A — > oo both ([TT| ) and ( p!4D reduce to 



Infill, (16) 



g2 



displaying the logarithmic divergence as A^ oo 



III. APPLICATION TO QUENCHED QCD 

The effects of quenching QCD may be categorized as perturbative or global effects and 
nonperturbative or channel specific effects. As discussed in the introduction, the effects in the 
perturbative regime may be accounted for through a simple renormalization of the coupling 
constant. However, the effects in the nonperturbative regime will, of course, be channel 
specific. The p-meson channel might be particularly vulnerable to the nonperturbative 
effects of quenching QCD due to the fact that the p is unstable in full QCD and stable 
in the quenched approximation. As we shall see however, our estimated "unquenching" 
corrections turn out to be rather small. 

Figure |^ displays p and squared pion masses for the three lightest quark masses used 
in the quenched QCD analysis of the GFll group for the 32 x 30 x 32 x 40 lattice at 
/3 = 6.17 for 219 configurations. The dashed line illustrates the linear relationship assumed 
in extrapolating the hadron masses to the critical point where the pion mass vanishes. 
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Provided the lattice spacing is defined by physical observables which do not have a similar 
dependence on pion decays, the physics associated with two-pion intermediate states of the 
p should simply be added onto the results extracted from calculations of quenched QCD. 
The correction is particular to the p meson and is unlikely to be accounted for if the lattice 
spacing is fixed by the nucleon mass for example. 

The solid line of figure | illustrates the addition of the self-energy correction to the linear 
extrapolation where the regulator mass = 1 GeV^ has been selected [|^. This choice of 
A is physically motivated and indicates the correction to the p mass at the physical point is 
negligible. The dot-dashed curve illustrates the correction when A^ = 2 GeV^. 

Figure ^ indicates that Ep^^ is less than 6% of the squared p mass for A^ < 2 GeV^. This 
corresponds to a less than 3% correction to the p mass itself. The conclusion that may be 
drawn from this analysis, which differs from previous considerations of these issues, is that 
the predictions of the GFll group are not "too good to be true". The channel specific 
nonperturbative correction to the p-meson mass may actually be rather small. 

The magnitude of the corrections estimated here is much smaller than that anticipated 
in the analysis of Ref. [^] for any reasonable choice of A. Geiger and Isgur predict an un- 
quenching correction to the p mass of approximately —160 MeV in contrast to our prediction 
of a to 25 MeV reduction of the p mass. 



IV. FULL QCD CALCULATIONS 

In full QCD simulations, the two-pion induced self-energy corrections are, of course, 
already included. However, the continuum predictions derived here will differ from those 



anticipated on the lattice |]T3[ largely due to the discretization of the momenta on the lattice 
and the lattice regularization itself. This renders the divergent integral of (|^) to a finite sum 
over a few two-pion states. In fact, if one hopes to recover the continuum physics, it will be 
necessary to first correct the hadron masses determined at unphysical quark masses to the 
continuum, infinite volume limit prior to extrapolating to physical quark masses. 
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Figure H suggests that it may be extremely difficult to see the effects of virtual two- 
pion intermediate states in full QCD. Since a great deal of the integral strength is lost for 
current lattice regularization parameters, the correction curve in figure ^ is most likely an 
optimistically large deviation from the linear relation. Even when full QCD calculations 
reach the current state of quenched calculations, it is questionable whether one will be 
capable of discerning the effects of virtual two-pion states of the p while the p is stable. Of 
course, once the p becomes unstable and decays to two pions, the lower lying two-pion states 
will need to be subtracted from the correlation function prior to extracting the p mass. This 
renders an examination of the p-meson mass above the two-pion cut nearly impossible. It 
should be mentioned that these arguments support the recent findings of Ref. [|17| where an 



attempt to observe the effects of virtual pion states in p correlation functions failed. 

Ultimately, full QCD calculations will reach the point where linearly extrapolating the p 
mass to physical quark masses ignoring nonanalytic effects will introduce a relevant system- 
atic error. To this end we present figure ^ which illustrates the relative amount which should 
be added onto the p mass extracted from a linear extrapolation of full QCD data corrected 
to the continuum, infinite volume limit |TB|. We remind the reader that these effects are 



being calculated using the dipole form for the cutoff which we have more or less arbitrarily 
chosen. Thus the curves are illustrative only. The x-axis indicates the point at which the 
derivative is determined for the linear extrapolation. Values for are as in figures |l| and |^. 
For example, if the effective physical value for is 1 GeV^ and the derivative is determined 
at ~ 0.25 GeV^, the p mass extracted from a linear extrapolation should be augmented 
by approximately 25 MeV prior to comparing with experiment. 

The size of the corrections to linear extrapolations illustrated in figure ^ are generally 
smaller than the 70 MeV addition predicted by Geiger and Isgur. However, it is possible 
to recover a correction to the p-meson mass extrapolation similar in magnitude to that of 
Geiger and Isgur's analysis, provided one linearly extrapolates from the onset of the cut as 
in their investigation. Of course, it is not practical to extrapolate from the onset of the cut 
in actual lattice calculations. In fact, figure ^ indicates that this is the worst possible place 
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to attempt an extrapolation to physical quark masses. The derivative is more likely to be 
averaged from a number quark masses corresponding to squared pion masses of 0.2 to 0.6 
GeV^. Hence the corrections are expected to be the order of 10 to 20 MeV and possibly 
neghgible if the effective physical value for is 2 GeV^ or more. 

An important point to mention here is that a simple calculation of the p-meson mass in 
full QCD will not circumvent the problems associated with its decay to pions. In fact, as 
one works hard to drive down the quark mass, linear extrapolations to the physical point 
will increasingly underestimate the p mass. Determination of the last few percent of the 
physical p-meson mass requires additional information describing the dependence of pm: 
interactions such that a suitable extrapolation function may be identified. Fortunately, these 
corrections are small and may be neglected until systematic uncertainties associated with 
the finite volume and finite spacings of the lattice are understood and eliminated. 

V. SUMMARY 

We have calculated the two-pion induced self-energy correction to the p-meson mass in 
a manner that allows an estimation of the correction to quenched QCD calculations and an 
analysis of extrapolations of full QCD results. 

The analysis of full QCD extrapolations indicates linear extrapolations of typical p-meson 
masses extracted from lattice correlation functions will underestimate the physical p mass 
by 10 to 20 MeV. An important point to draw from the analysis is that as the quark masses 
become lighter, linear extrapolations to the physical point will increasingly underestimate 
the p-meson mass. 

We estimate the corrections to quenched calculations to be the order of a few percent 
and quite possibly neghgible (0 to —25 MeV). These results lend credence to the success of 
quenched QCD in describing the physical low-lying hadron mass spectrum. 
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FIGURES 

FIG. 1. The real part of SpTTTr from ([Tl| ) evaluated at = m?p for a variety of cutoffs, at 
s = A^. In this and the following figures, the finely dashed vertical line marks the position of the 
physical point. For small cutoffs, most of the strength in the integral lies below the p mass and as 
a result the p mass is pushed up. 

FIG. 2. The two-pion induced self-energy (a) and its derivative with respect to (b) for 
a regulator mass of = 2 GeV^. Both real (solid line) and imaginary (dashed line) parts are 
illustrated. The derivative clearly displays the nonanalytic behavior encountered at rup = 2mj^. 

FIG. 3. The real part of the self-energy for dipole regulator masses A taking the same values 
used in figure ||. The sensitivity of the results to the value of A is greatly reduced. 

FIG. 4. p and squared pion masses for the three lightest quark masses used in the quenched 
QCD analysis of the GFll group on their largest lattice. The dashed line illustrates the linear 
relationship assumed in extrapolating the hadron masses to the critical point. The solid curve 
displays the self-energy correction for a regulator mass of A^ = 1 GeV^ where the correction to 
the p mass at the physical point is negligible. The dot-dash curve corresponds to A^ = 2 GeV^. A 
lattice scale parameter of 2.73 GeV has been applied to the otherwise dimensionless lattice results. 

FIG. 5. The amount to be added to the p mass extracted from a linear extrapolation of full 
QCD data. The x-axis indicates the point at which the derivative is determined for the linear 
extrapolation. Values for A^ are as in figures |^ and ^. 
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at ~ 0.25 GeV^, the p mass extracted from a linear extrapolation should be augmented 
by approximately 25 MeV prior to comparing with experiment. 

The size of the corrections to linear extrapolations illustrated in figure 5 are generally 
smaller than the 70 MeV addition predicted by Geiger and Isgur. However, it is possible 
to recover a correction to the /j-meson mass extrapolation similar in magnitude to that of 
Geiger and Isgur's analysis, provided one linearly extrapolates from the onset of the cut as 
in their investigation. Of course, it is not practical to extrapolate from the onset of the cut 
in actual lattice calculations. In fact, figure 5 indicates that this is the worst possible place 
to attempt an extrapolation to physical quark masses. The derivative is more likely to be 
averaged from a number quark masses corresponding to squared pion masses of 0.2 to 0.6 
GeV^. Hence the corrections are expected to be the order of 10 to 20 MeV and possibly 
negligible if the effective physical value for is 2 GeV^ or more. 

An important point to mention here is that a simple calculation of the /j-meson mass in 
full QCD will not circumvent the problems associated with its decay to pions. In fact, as 
one works hard to drive down the quark mass, linear extrapolations to the physical point 
will increasingly underestimate the p mass. Determination of the last few percent of the 
physical /j-meson mass requires additional information describing the dependence of /jtttt 
interactions such that a suitable extrapolation function may be identified. Fortunately, these 
corrections are small and may be neglected until systematic uncertainties associated with 
the finite volume and finite spacings of the lattice are understood and eliminated. 

V. SUMMARY 

We have calculated the two-pion induced self-energy correction to the /j-meson mass in 
a manner that allows an estimation of the correction to quenched QCD calculations and an 
analysis of extrapolations of full QCD results. 

The analysis of full QCD extrapolations indicates linear extrapolations of typical /j-meson 
masses extracted from lattice correlation functions will underestimate the physical p mass 
by 10 to 20 MeV. An important point to draw from the analysis is that as the quark masses 
become lighter, linear extrapolations to the physical point will increasingly underestimate 
the /j-meson mass. 

We estimate the corrections to quenched calculations to be the order of a few percent 
and quite possibly negligible (0 to —25 MeV). These results lend credence to the success of 
quenched QCD in describing the physical low-lying hadron mass spectrum. 
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continuum, infinite volume limit prior to extrapolating to physical quark masses. 

Figure 4 suggests that it may be extremely difficult to see the effects of virtual two- 
pion intermediate states in full QCD. Since a great deal of the integral strength is lost for 
current lattice regularization parameters, the correction curve in figure 4 is most likely an 
optimistically large deviation from the linear relation. Even when full QCD calculations 
reach the current state of quenched calculations, it is questionable whether one will be 
capable of discerning the effects of virtual two-pion states of the p while the p is stable. Of 
course, once the p becomes unstable and decays to two pions, the lower lying two-pion states 
will need to be subtracted from the correlation function prior to extracting the p mass. This 
renders an examination of the /j-meson mass above the two-pion cut nearly impossible. It 
should be mentioned that these arguments support the recent findings of Ref. [17] where an 
attempt to observe the effects of virtual pion states in p correlation functions failed. 

Ultimately, full QCD calculations will reach the point where linearly extrapolating the p 
mass to physical quark masses ignoring nonanalytic effects will introduce a relevant system- 
atic error. To this end we present figure 5 which illustrates the relative amount which should 
be added onto the p mass extracted from a linear extrapolation of full QCD data corrected 
to the continuum, infinite volume limit [18]. We remind the reader that these effects are 
being calculated using the dipole form for the cutoff which we have more or less arbitrarily 
chosen. Thus the curves are illustrative only. The x-axis indicates the point at which the 
derivative is determined for the linear extrapolation. Values for are as in figures 1 and 3. 
For example, if the effective physical value for is 1 GeV^ and the derivative is determined 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 

m/ (GeV^) 

FIG. 5. The amount to be added to the p mass extracted from a linear extrapolation of full 
QCD data. The s-axis indicates the point at which the derivative is determined for the linear 
extrapolation. Values for A^ are as in figures 1 and 3. 
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FIG. 4. and squared pion masses for the three lightest quark masses used in the quenched 
QCD analysis of the GFll group [3] on their largest lattice. The dashed line illustrates the linear 
relationship assumed in extrapolating the hadron masses to the critical point. The solid curve 
displays the self-energy correction for a regulator mass of = 1 GeV^ where the correction to 
the p mass at the physical point is negligible. The dot-dash curve corresponds to = 2 GeV^. A 
lattice scale parameter of 2.73 GeV has been applied to the otherwise dimensionless lattice results. 

Figure 3 indicates that Sp^^ is less than 6% of the squared p mass for A^ < 2 GeV^. This 
corresponds to a less than 3% correction to the p mass itself. The conclusion that may be 
drawn from this analysis, which differs from previous considerations of these issues, is that 
the predictions of the GFll group [3] are not "too good to be true". The channel specific 
nonperturbative correction to the /j-meson mass may actually be rather small. 

The magnitude of the corrections estimated here is much smaller than that anticipated 
in the analysis of Ref. [5] for any reasonable choice of A. Geiger and Isgur predict an un- 
quenching correction to the p mass of approximately —160 MeV in contrast to our prediction 
of a to 25 MeV reduction of the p mass. 



In full QCD simulations, the two-pion induced self-energy corrections are, of course, 
already included. However, the continuum predictions derived here will differ from those 
anticipated on the lattice [13] largely due to the discretization of the momenta on the lattice 
and the lattice regularization itself. This renders the divergent integral of (7) to a finite sum 
over a few two-pion states. In fact, if one hopes to recover the continuum physics, it will be 
necessary to first correct the hadron masses determined at unphysical quark masses to the 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

FIG. 3. The real part of the self-energy for dipole regulator masses A taking the same values 
used in figure 1. The sensitivity of the results to the value of A is greatly reduced. 



III. APPLICATION TO QUENCHED QCD 



The effects of quenching QCD may be categorized as perturbative or global effects and 
nonperturbative or channel specific effects. As discussed in the introduction, the effects in the 
perturbative regime may be accounted for through a simple renormalization of the coupling 
constant. However, the effects in the nonperturbative regime will, of course, be channel 
specific. The /j-meson channel might be particularly vulnerable to the nonperturbative 
effects of quenching QCD due to the fact that the p is unstable in full QCD and stable 
in the quenched approximation. As we shall see however, our estimated "unquenching" 
corrections turn out to be rather small. 

Figure 4 displays p and squared pion masses for the three lightest quark masses used 
in the quenched QCD analysis of the GFll group [3] for the 32 X 30 X 32 X 40 lattice at 
/3 = 6.17 for 219 configurations. The dashed line illustrates the linear relationship assumed 
in extrapolating the hadron masses to the critical point where the pion mass vanishes. 

Provided the lattice spacing is defined by physical observables which do not have a similar 
dependence on pion decays, the physics associated with two-pion intermediate states of the 
p should simply be added onto the results extracted from calculations of quenched QCD. 
The correction is particular to the p meson and is unlikely to be accounted for if the lattice 
spacing is fixed by the nucleon mass for example. 

The solid line of figure 4 illustrates the addition of the self-energy correction to the linear 
extrapolation where the regulator mass A^ = 1 GeV^ has been selected [16]. This choice of 
A is physically motivated and indicates the correction to the p mass at the physical point is 
negligible. The dot-dashed curve illustrates the correction when A^ = 2 GeV^. 
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FIG. 2. The two-pion induced self-energy (a) and its derivative with respect to (b) for 
a regulator mass of = 2 GeV^. Both real (solid line) and imaginary (dashed line) parts are 
illustrated. The derivative clearly displays the nonanalytic behavior encountered at mp = 2m^. 
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This form maintains the normahzation of (/p^^ defined at the physical p mass and renders 
the integral of (8) finite. 

Of course, this approach is not without a few unphysical side effects. The most obvious 
problem is the introduction of spurious poles in the space like region for the s dependence 
of (/pjTTr- However, the dispersion integral only samples the time like region and the presence 
of these unphysical poles should not affect the results. One could consider other functional 
forms. However, the effective physical value for the regulator mass. A, is itself unknown. 
Our aim is to estimate the importance of the two-pion induced self-energy relative to the p 
mass, as opposed to attempting to evaluate the actual correction. For this reason we view 
a consideration of the dipole regulator to be adequate. 

Evaluation of the dispersion relation of (8) with (13) leads to the following functional 
form for the real part of Sp^^ 

8m2 12m2 

lOm^ 6m2 o2\ ( R^-l\ 

+ |l + ^ + ^Uln 1^ (14) 



where 




Re S„. = -;^ ! + + 





^ 2 \ 1/2 / ^ 2 \ 1/2 



J , and /3a= l^l + ^j . (15) 

The imaginary part is recovered as in (9). 

Figure 2 illustrates the real part of the self-energy and its derivative with respect to 
at = 2 GeV^. The derivative clearly displays the nonanalytic behavior encountered at 
nip = 2m^. The second derivative is discontinuous at nip = 2m^ and is infinite from above. 
The imaginary part of the self-energy is also illustrated in figure 2. 

A comparison with figure 1 indicates that at A^ = 2 GeV^ the results are not too sensitive 
to the manner in which the integral is regulated. Figure 3 illustrates the real part of the 
self-energy for the same values of A^ used in figure 1. The sensitivity of the results to the 
value of A is greatly reduced and all curves display the same qualitative behavior. 

In the limit of A — )• oo both (11) and (14) reduce to 

4:0 /1- I (/- \ H' J ' \^ ~^ q / 

(16) 

displaying the logarithmic divergence as A^ — )• oo. 
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FIG. 1. The real part of SpTrvr from (11) evaluated at = for a variety of cutoffs, at s = A^. 
In this and the following figures, the finely dashed vertical line marks the position of the physical 
point. For small cutoffs, most of the strength in the integral lies below the p mass and as a result 
the p mass is pushed up. 
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Figure 1 illustrates the real part of Sp^^ evaluated at 



m. 



for a variety of cutoffs 

ranging from slightly above to 4 GeV^. For small cutoffs, most of the strength in the 
integral lies below the p mass and consequently the p mass is pushed up due to the mixing 
with pion states. Of course, this behavior is completely consistent with that anticipated by 
simple quantum mechanical arguments. For larger cutoffs the strength above the p mass 
acts to reduce the p mass. 



B. Dipole Cutoff 

While the ^-function is useful as an illustrative tool, it suffers from being physically 
artificial and the results can be very sensitive to the value of the cutoff, as illustrated in 
figure 1. In an attempt to better represent the q^ dependence of gp^y^y a monopole form for 
each vertex is introduced and the dispersion relation of (8) is evaluated with 

2 (l±JlV 

A'pTTTT 9p7T7T I ^ _^ ^2 I " ^^^f 
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^ S,.. [g^^ - ^ j . (6) 
Y^TT -g gjven by the standard one loop integrals 

_-v-= /J!^ 2 [ (g--2r)(g--2r) 2g,. ] 

y (2 7r)^ ^"^^ [ [(^ - kf -ml + ze] [F - ml + ze] k^-ml + zej ' ^ > 



Physically, the integral is convergent due to the momentum dependence of (/pTrvr- However, 
this momentum dependence is unknown. In light of this uncertainty, it is reasonable to 
parameterize the momentum dependence in terms of some given functional form with an 
adjustable parameter controlling how the function falls off as a function of momentum trans- 
fer. We consider two regulation prescriptions. In one, we assume a monopole form for (/pTrvr, 
and for comparison, we also consider a sharp ^-function cutoff. 

The simplest fashion for introducing a covariant cutoff function is through the use of a 
dispersion relation. The second term of (7) is q independent and serves only to subtract the 
quadratic divergence of the first term in maintaining current conservation. As a result, we 
write a dispersion relation for ^{q'^) with one subtraction at q^ = 0, 

W)^- ds^ 8 

TT Jo s s — q'^ 

Of course, the imaginary part of Sp^^ may be easily determined using any number of tech- 
niques for rendering the integral of (7) finite. The imaginary part is 

ImSp..(,^) = |^,^ {l-^-^)''' 0[q'-^ral) . (9) 

The value of (/p^^ at q^ = is fixed by equating the imaginary parts of 

+ Ep^^ = (^mp + , (10) 

at q^ = . The physical values [15] nip = 768.1 MeV and V = 151. 5 MeV fix (/p^^ at ~ 6.0. 



A. ^-Function Cutoff 

To illustrate the physics associated with the real part of the self-energy we first consider 
the integral of (8) cut off covariantly by a sharp ^-function at s = A^. The functional form 
is 

= 1^ e (in (i-^) + ^ .^ - In (^) | . (11) 

48 TT^ [ W + at^J q^ ycfk + CFq) J 

where 
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bosons such as the pion, as they do not incorporate chiral symmetry. While there are models 
of the piTTT vertex which incorporate chiral symmetry and chiral symmetry breaking [10-12], 
these models are based on particular dynamical assumptions. Accordingly, it is difficult to 
assess the reliability of such models. Instead, we consider a range of possibilities for the 
vertex and to simplify this task, we consider convenient phenomenological forms. 

One other paper addressing this issue [13] sidesteps the problems surrounding the de- 
pendence of QpTy^y by fixing (/p^^ to a constant and making two subtractions of the divergent 
integral at (j'^ = 0. These subtractions are absorbed into a mass and wave function renormal- 
ization. However, this approach excludes any analysis of /j-meson mass extrapolations as the 
subtraction terms themselves have an unknown dependence which has been lost in the 
renormalization procedure. Moreover, contributions from virtual two-pion states have been 
absorbed into the bare lattice parameters which is inconsistent with the dynamics contained 
in the quenched approximation. 

The outline of this paper is as follows. In Section II the model used in examining the 
two-pion induced self-energy is outlined. Two methods for regulating the divergent self- 
energy are explored. In section III the relevance of the self-energy corrections to quenched 
QCD simulations is discussed. Section IV addresses the quark mass extrapolation of full 
QCD calculations and the importance of nonlinear behavior in the /j-meson mass. Finally, 
the implications of this investigation are summarized in Section V. 



II. THE SELF-ENERGY 



In modeling the two-pion induced self-energy of the p meson, Sp^^, the standard /jtttt 
interaction motivated by low-energy current algebra is used. The effective Lagrange inter- 
action has the form [14] 



9p^^ P" (tT TT^ + gj^^ TT^ p"^ . 



The pions are further assumed to interact exclusively through the p channel as summarized 
in the following Schwinger-Dyson equation for the p propagator 



G/iiy — 6*°^, -|- (7°^ E*^^ Griy , (3) 



where 



^A^- - q2_ + le y^'"' J ' 

in Landau gauge, and Mq is the bare /j-meson mass. The self-energy Sp^^ is defined through 
the solution 

^ - Mi - Ep.. + ze {^'^ ~ ' 

where 
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self-energy of the p meson in relation to lattice QCD calculations [5]. Their results are based 
on a string breaking quark model and predict large corrections to smooth extrapolations of 
the /j-meson mass at approximately 70 MeV. These authors were motivated by the long 
standing problem of QCD predictions for the A^//j-mass ratio being too large. Their hope 
was the two-pion induced self-energy correction would sufficiently raise the p mass to solve 
this problem. However, we now understand that both the finite lattice spacing and the finite 
volume of the lattice act together to push up this mass ratio. Current estimates [3] for this 
ratio corrected to the infinite volume, continuum limit are f .22ib0.f f in excellent agreement 
with the experimental value of f .22. 

Geiger and Isgur [5] advocate using the nonlinearities in the p mass as a function of the 
quark mass to correct for the linear extrapolation of lattice results. We wish to stress that 
such a procedure is only sensible for the extrapolation of full QCD calculations. As we shall 
argue, the entire two-pion induced self-energy is absent in the quenched approximation and 
should be added on to quenched QCD results prior to comparison with experiment. We 
note that for the model of Ref. [5], quenched lattice calculations of the /j-meson mass would 
be reduced by 160 MeV instead of increased by 70 MeV. This would further exacerbate the 
^^m]\j jnip problem" discussed in their paper, rather that curing it. 

In the quenched approximation, the p meson cannot decay to two light pseudoscalars. 
As discussed in Ref. [6] it is not possible to generate intermediate states of the p meson in 
the quenched approximation in which one has two isovector pseudoscalars. One might worry 
about the presence of the light isoscalar pseudoscalar rj' which fails to obtain its heavy mass 
in the quenched approximation [6-9]. However, decay of the p meson to a tt ?y' is forbidden 
by G-parity and decay to two ?y'-mesons is of course forbidden by the isospin invariance of 
the strong interactions. Hence the physics associated with two-pion intermediate states of 
the p should simply be added onto the results extracted from calculations of quenched QCD, 
provided the lattice spacing is defined by physical observables which do not have a similar 
dependence on pion decays. 

Current calculations of full QCD typically employ quark masses which place 2m^ > 
nip. As a result, the functional form of the extrapolation function should account for the 
nonanalytic behavior in the /j-meson mass as the two-pion cut is encountered. On the lattice, 
the spectral density does not have a cut but rather a series of poles at the points satisfying 



where pn are the discrete momenta allowed on the lattice. Obviously, to fully account for 
the two-pion induced self-energy of the /?, one must first extrapolate the lattice results to 
zero lattice spacing and infinite volume prior to extrapolating the quark masses to physical 
values. 

In evaluating the integrals describing the coupling of pions to the p one must take into 
account the dependence of the p to two tt coupling constant (/p^^ refiecting the internal 
structure of these mesons. While this dependence was extracted from a string breaking 
quark model in Ref. [5], we elect to take a more agnostic approach and consider different 
methods of cutting off the integral. In particular we investigate a sharp ^-function cutoff and 
a dipole cutoff. The underlying reason for our agnosticism is the belief that models of the 
sort underlying Ref. [5] are not likely to correctly describe the structure of pseudo-Goldstone 
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I. INTRODUCTION 



The lattice regularized approach to quantum field theory provides the best forum for the 
examination of the fundamental nonperturbative aspects of QCD. In the low momentum 
transfer regime, it is the only approach which in the foreseeable future holds a reasonable 
promise of confirming or rejecting the validity of QCD as the underlying theory of the strong 
interactions. 

Recently, it has become possible to perform quenched QCD calculations in which all 
systematic uncertainties are quantitatively estimated. Thus if the effects of quenching can 
be understood, the validity of QCD may be tested in the nonperturbative regime. Of 
particular note is the recent determination of the QCD coupling constant, aj^^ from the 
IS — IP mass splitting of charmonium [I]. In this case one believes that the effects of 
quenching may be estimated with minimal model dependence. Corrections to the utilization 
of the quenched approximation of QCD have been estimated and are currently the dominant 
source of uncertainty in the final predictions [2]. 

In this paper we will continue efforts along this line through the examination of systematic 
uncertainties in hadron mass spectrum calculations. In particular, the importance of the 
two-pion induced self-energy contribution to the /j-meson mass is evaluated in relation to 
the quenched approximation of QCD and to full QCD. In the quenched approximation, the 
physics associated with two-pion intermediate states has been excluded in the numerical 
simulations. 

This investigation is motivated by recent results from the GFII group [3] for the low-lying 
hadron mass spectrum in the quenched approximation of QCD. Their analysis is the first 
to systematically extrapolate QCD calculations to physical quark mass, zero lattice spacing 
and infinite volume. Their predictions display an impressive agreement with experiment. 
Of eight hadron mass ratios, six agree within one standard deviation and the remaining two 
ratios agree within 1. 60". However, since the quenched approximation leaves out so much 
important physics, one might question whether these results are actually too good [4]. 

In the perturbative regime, many of the effects of not including disconnected quark loops 
when preparing an ensemble of gauge configurations may be accounted for in a simple renor- 
malization of the strong coupling constant. However, one also anticipates nonperturbative 
effects in making the quenched approximation. Unlike a global renormalization of the cou- 
pling constant, these effects are expected to be channel specific. For example, the quenched 
approximation of QCD leaves out the physics associated with the decay of the p meson to 
two pions. This physics must be accounted for and added to the quenched results prior to 
comparing with experimental data. Moreover, the calculated hadron masses are extrapo- 
lated as a function of the pion mass squared to the point at which the pion mass vanishes 
using linear extrapolation functions. Such an approach neglects nonlinear and indeed, non- 
analytic behavior in the continuum extrapolation function. For example, in the case of the 
/j-meson mass, one expects nonanalyticity associated with the onset of the two-pion cut. 

A priori one does not know the relative importance of two-pion intermediate states of 
the p meson in describing the p mass. The substantial width of the p meson at 151. 5 it 1.2 
MeV indicates its coupling to pions is not small and correspondingly these dynamics may 
have significant infiuence on the /j-meson mass. 

Geiger and Isgur were the first to study the possible importance of the two-pion induced 
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Abstract 



Two-pion induced self-energy contributions to the />-meson mass are exam- 
ined in relation to the quenched approximation of QCD, where the physics 
associated with two-pion intermediate states has been excluded from vector- 
isovector correlation functions. Corrections to quenched QCD calculations of 
the />-meson mass are estimated to be small at the order of a few percent of 
the />-meson mass. The two-pion contributions display nonanalytic behavior 
as a function of the pion mass as the two-pion cut is encountered. The impli- 
cations of this nonanalytic behavior in extrapolations of full QCD calculations 
are also discussed. We note that for full QCD, the error made in making a 
linear extrapolation of the p mass, neglecting nonanalytic behavior, increases 
as one approaches the two-pion cut. 



